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PERTURBATION- ITERATION THEORY 


FOR ANALYSING MICROWAVE STRIPLINES 

ABSTRACT 

by 

BRIAN EDWARD KRETCH 

A perturbation-iteration technique is presented for determin- 
ing the propagation constant and characteristic impedance of an 
unshielded microstrip transmission line. The method converges to 
the correct solution with a few iterations at each frequency and 
is equivalent to a full wave analysis. The perturbation-iteration 
method gives a direct solution for the propagation constant with- 
out having to find the roots of a transcendental dispersion 
equation. 

The theory is presented in detail along with numerical 
results for the effective dielectric constant and characteristic 
impedance for a wide range of substr^^e dielectric constants, 
stripline dimensions, and frequencies. 
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CHAPTER I 


INTRODUCTION 

1.1 Introduction 

The analysis of Infinite, straight, unshielded microwave strip 
lines Is an area of great Importance In microwave field theory. 
Microwave striplines are used extensively in high frequency communi 
cation networks fcr radar and satellite systems. Figure 1.1 shows 
the cross section of typical unshielded and shielded microstrip 
lines and also the related slot line. The effectiveness of 
their use in circuits is contingent on how well the designer can 
understand the striplines to know their properties and character- 
istics. There has been a great deal of previous research in the 
area of modeling and analysis of microwave striplines that encom- 
passes three basic approaches: 1) static or quasi-TEM analysis, 

2) dispersion modeling, and 3) fullwave analysis. 

Wave propagation along a microwave stripline is similar to 
pure TEM propagation. The differences arise due to the presence 
of the ground plane and the dielectric substrate. At low frequen- 
cies, these deviations are negligable so that the waves can be 
analyzed as pure TEM waves. This is the foundation of the static 
or quasi-TEM analysis. As the frequency increases, the frequency 
effects become more predominant and the analysis must also include 
these effects. This is the basis of the fullwave analysis. Models 
of the microwave stripline can also be developed to simplify the 
evaluation of high frequency properties. Such is the method used 
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Figure 1.1 (a) Unshielded Microstrip Line, (b) Shielded Microstrip 

Line, (c) Slot Line, (d) Shielded Slot Line. 
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in dispersion models. 

The different approaches for analyzing the microwave stripline 
all utilize the same basic boundary conditions. These boundary 
conditions all emanate from Maxwell's Equations; 

7 x f ■ -jw$ 

7 x t - jwS + 5 
7 • B • p 

7*1-0 Cl- 1) 


The media will be assumed to be isotropic which permits the following 
relations : 

D - ice ^ 
o 

5 - ji 5 (1.2) 

o 

The boundary conditions will thus encompass the following; 



at y - 0 and on the mlcrostrlp 

are continuous across the air-dielectric 
interface 


(1.3) 


where and are the components of the total current density 

and p is the total charge density on the microstrip. 
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In Sec. 1.2 - 1.4, a brief review of past work on microwave 
stripline propagation will be presented. This review will be 
brief with further references provided so that a more complete 
understanding of them may be acquired If desired. 

1. 2 Quasi- T.H Analysis 

The first type of analysis examined will be the static, or 
quasi-TEM analysis. As previously mentioned, the wave propagation 
at low frequencies (up to 1GHz) can be approximated by a pure TEM 
wave. The foundation of this analysis Is based on the previously 
stated boundary conditions and contingent on a few basic assumptions. 
Since the frequency, u , Is assumed to be small and approach zero, 
the terms jojff and juilt iu c’axwell's equations (1-2) may be 
assumed to equal zero. In a different sense, this assumption states 
that as a) approaches zero, the wavelength X approaches infinity. 
This assumption introduces negligible error into the stripline 
analysis at sufficiently low frequencies. Thus, at any location 
along the microstrip but not on the strip, 

7 x £ - -juy 5^0 

O 

7 x 5 - JwS % 0 

On the strip, the condition still demands that 

7x5- 


Using the relationship 


5 


t - 7 x t 

along with Che Lorentz condition 

7«J « “j«u e >c$ 
o o 

Che following two relationships can be developed: 
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These assumptions and equations form the basis by which different 
quasi-static analyses are developed and employed. 

The essential quantities to be solved for in the quasi- TEM 
analysis are the characteristic impedance Z fi and the propagation 
constant 8 . These entities are obtained by first determining two 
static capacitances p?r unit length along the stripline. The first 
capacitance C is found by replacing the substrate by air while 
maintaining the same configuration. The second capacitance C q 
is defined as that of the stripline with the substrate present. 

With these two quantities, the values of Z ; and 8 may be solved 

for by using the following equations: 


Z - (1/c)/ 1/C C 
a o 


8 


k 

o 


/ C /C 
o a 


(1.4a) 


(1.4b) 
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where 

c - Speed o£ light 



An alternative approach to solving for and g Is to de- 

dermine the static Inductance per unit length of the mlcrostrlp 
L q using the relationship 

L I « A on the microstrip 
o o z r 

where 

I ■ total z-directed current 

o 

A ■ z-dlrected vector potential 

and leading to the solutions 

8 -a) / L C 

o o 



From the solution for the propagation constant 8 , an effective 
dielectric constant k may be obtained. This is given by 


'eff 


fcf ) 2 


The significance of is to have a simple direct relationship 

between the quasi-TEM wave with a dielectric and the quasi-TEM wave 
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with a dielectric having an effective dielectric constant *< e f£* 

This permits easy calculations of 6 at various frequencies. 

There are many techniques used for solving for and 8 
along the infinite line. The following summarizes the more prominent 
methods and lists references for further research if desired. 

The first method is the Conformal Transformation Method [1-3] 

This method solves for C and C using a conformal trans forma- 

o a 

tlon. The difficulty with this method is that the viicrostrlp 
structure cannot be conformally mapped into a different fom which 
is easily analyzed. Approximations would have to Le incorporated 
into the analysis which would decrease the accuracy of the method. 

The second type of static analysis is the Variational Method [4-5] 

Ir this method, a stationary expression for the capacitance in terms 
of charge distribution is formed. The charge distribution is 
represented by a functional fom containing one or more variational 
parameters. The best approximation to the capacitance is then ob- 
tained by choosing the variational parameters so as to make the 
expression for the capacitance stationary, for a given accuracy, 
the computations are of the same order of complexity as that of 
other methods. The third method of analysis is the Relaxation 
Method [6 ]• This technique divides the cross section of the strip- 
line into a grid, assuming values of potentials at all points, 
and then modifying these values, or "relaxing" them, by using a 
finite difference approximation to Laplace’s equation for the 
potential field. The final method is termed the Integral Equation 
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Method t 7 ]. This technique uses a Green's Function to formulate 
an integral equation to solve for the charge density. The capaci- 
tances are then determined from the charge density and potential, 
which has been presumed to have assigned a constant value on the 
conducting strip. These methods are the main ones employed in 
solving for the unknown Z £ and 8 under static conditions. 

1. 3 Dispersion Models 

The second type of analysis of the straight. Infinite, 
unshielded microwave stripline is developing and applying disper- 
sion models. The static analysis is useful in that it is an easy 
approach for studying the waves at low frequencies, Its principal 
deficiency, though, results from its failure to account for the 
non-TEM nature of the waves at higher frequencies. In particular, 

8 is not a linear function of ui so the stripline is dispersive. 
Also, the characteristic Impedance increases slowly with frequency. 
At sufficiently high frequencies, other modes begin to propagate. 
Consequently, various ad hoc dispersion models have been developed 
which help to analyze the striplines and take into account these 
non-TEM characteristics. Furthermore, they help to present direct 
relationships for and in accordance with a specific 

frequency. These models can thus contribute to a better understand- 
ing of the wave propagation along the stripline over a more diverse 
range of frequencies. 

The first type of dispersion model is the Coupled TEM and TM 
modes Model [ 8 ]. This model developes a direct equation for 
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based on the assumption that the non-TJSM characteristics can be 
attributed to the lowest order surface wave TM q mode. Empirical 
Relations [ 9 ] have also been developed, based on various theor- 
etical and experimental studies of wave dispersion which yield a 
direct equation for the wave phase velocity. The Dielectric- 
loaded Ridged Waveguide Model [10] and the Planar Waveguide Model 
[11] alter the configuration of the microwave stripline structure 

for ease of direct mathematical computation of ic and Z while 

eft c 

maintaining the same dispersion characteristics as that of the 

strlpllne. The final dispersion model is the Coupled Transmias on 

Model [12], T.»is model describes the stripline in terms of coupled 

TEM and TE modes using coupled transmission lines and basic circuit 

analysis to solve for K e f£ and Z . The consequences of using 

dispersion models will be to cause the strlpllne to be easier to 

mathematically analyze while sacrificing some accuracy in the results. 

The various approximate dispersion models are useful in many 

engineering applications but generally are not sufficiently accurate 

in critical applications where the exact values of Z and $ are 

c 

needed. 

1.4 Full Wave Analysis 

The final and most accurate microstrip analysis is the Fullwave 

Analysis [13-17]. This method makes no quasi-static assumptions while 
analyzing the microstrip configuration exactly as given with no con- 
version to analytically simpler, but less accurate, dispersion models. 
The fjllwave analysis employs an integral equation method in either 
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the space or Fourier transform domain and Is solvec using the method 
of moments • These different techniques are Initially based on equiva- 
lent representations of the waves along the microstrip. The analysis 
then assumes that the propagating wave is composed of hybrid TE and 
TM modes along the line due to two separate scalar potentials. The 
individual wave entities due to separate modes are converted to the 
Fourier Transform Domain (FTD) where they can be evaluated. By 
applying the correct boundary conditions In conjunction with one of 
the af foment lone d methods of analysis, the waves can be solved for 
along a stripline. The houndary conditions used are those given in 
Sec. 1.1 with no assumptions made concerning the waves or strip- 
line. The spatial domain solution applies the boundary conditions 
using a Green's function integral in the space domain. The unknown 
current is expanded in terms of a suitably chosen set of basis 
functions. This will result in computation of the amplitudes of the 
longitudinal and transverse currents using the method of moments. 

The FTD solution uses a Green' s function applied in the FTD to solve 
for the current amplitudes. This is done using Galerkin's Method 
along with an infinite weighted set of current density basis func- 
tions. The fullwave analysis can thus be used to examine the wave 
nature of the stripline with extreme accuracy with one important 
consequence: the analysis is very complicated with a great deal of 

n ume rical analysis required. The propagation constant 3 comes from 
the dispersion equation that is obtained and since this dispersion 
equation is a complicated transcendental equation, finding its roots 


a <- 


n 


is quite complex. The dispersion equation gives the propagation 
constant of higher inodes as well as that of the quasi-TEM mode. 

This important factor must be taken into account when choosing a 
technique for otripline analysis. 

1.5 Perturbation- Iteration Method 

The methods for analyzing the infinite straight, unshielded 
microwave stripline that have been previously developed each have 
their useful lness as well as their deficiences. The quasi-TEM 
analysis is adequate for analyzing the microstrip at low frequen- 
cies in a simple manner but is not sufficiently accurate at high 
frequencies. The dispersion model permits the analysis to be done 
with mathematical ease while sacrificing accuracy due to the 
approximations inherent to the models. The fullwave analysis is 
very accurate at all ranges of frequencier but is extremely cumber- 
some in computations. A new technique has been developed to 
analyze the waves along the stripline which offers greater simpli- 
city than existing fullwave analysis. This method, the Perturbation- 
Iteration Method, initially solves for the field components at quasi- 
TEM frequency values. As the frequency increases, the entitles are 
found using frequency dependent Green* s functions and an iteration 
method to successively yield greater accuracy. At each iteration 
level the value of 8 as the previously determined one is used to 
formulate the new frequency dependent Green's functions at the new 
frequency. The Iteration is then carried out to yield a new more 
accurate 6 and a field solution that satisfies the field equations 


intXa *4vM U 
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at higher frequencies. The effective frequency range is theoret- 
ically infinite which makes this technique attractive and useful. 

The most important feature that distinguishes this method from the 
full wave analysis described above is that £ is found explicitly 
instead of being determined from a complex dispersion equation. 

The pertinent equations and techniques of the Perturbation-Iteration 
Method will be developed for the complete analysis of the waves 
along the strip line in this thesis. The theory is developed in 
Chap. 2 along with the required Green's functions. Chapter 3 
presents extensive numerical and graphical results for the propa- 
gation constant & as a function of dielectric constant of the 
substrate, width to height ratio of the stripline, and frequency. 

The corresponding results for characteristic impedance are also 
presented. The final chapter is the concluding one and contains 
suggestions for further work and a brief discussion of a pertur- 
bation method that was initially developed but abandoned in favor 
of the Perturbation-Iteration Method. 


CHAPTER 2 


PERTURBATION-ITERATION METHOD THEORY 
2. 1 Introduction 

The Perturbation-Iteration Method used for analyzing the 
Infinite, straight, unshielded, microwave stripline is a 
straightforward approach which utilizes two static Green's 
functions, three frequency dependent Green's functions, and 
Fourier representations of the quantities associated with 
the raicrostrip structure. As the theory is developed, certain 
assumptions and approximations are made to help carry out the 
procedure in a much easier and faster far.hlon. The stripline 
structure is displayed again in Pig. 2.1 with an important 
assumption applied to it: ground plane sidewalls are put on 

both sides of the microstrip structure to allow for Fourier 
representation of the fields along the stripline. The valid- 
ity of inserting the side walls is based on the assumption 
that all of the fields decay to approximately zero at a reason- 
able distance from the microstrip. The stripline is oriented 
to allow for propagation of the fields in the +z direction 
with a propagation constant 0 . The dielectric substrate is 
considered to be uniform, linear, isotropic and possess a 
dielectric constant k . The thickness of the stripline is 
taken to be infinitely thin and therefore neglected. For 
certain calculations, the stripline i3 assumed to be just 
above the dielectric substrate in the air-filled region at a 
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Figure 2.1 Unshielded Microstrip Line with Ground Plane Side Walls. 
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height of y “ H + . 

The basic approach of the Perturbation- Iteration Method 

begins with a static analysis of the fields by letting u • 0. 

This yields a static charge density p and +2 directed 

current density J being determined and used to solve for 8 

z 

and the characteristic Impedance Z . The next step is to 

c 

Increase the frequency by a small Increment allowing 8 to 

be approximated from the static analysis. New values of p 

and J z are subsequently determined which are used to solve 

for new 8 and Z at the given frequency and which are 

c 

more accurate than the previously approximate values. This 
scheme may be repeated while maintaining u> constant which 
improves the values of 8 and Z or by increasing <a to 
yield new values of 8 and Z^ so as to cover a broader 
range of frequencies. At any one frequency, three or four 
iterations are usually sufficient to give values for 6 and 
Z c that are essentially the exact converged values. 

Before the mathematical theory is developed, an enumer- 
ation procedure is established to clarify the interpretation 
of the various field quantities. This theory involves a 
continuously iterating process which must be monitored so 
that proper field entity values of particular iterative 
levels are operated on correctly. The various field entities, 
such as E, H, 1?, $, and p, all contain a subscript to 

signify the level of iteration being referred to. As an 
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example, 4 has an Iterating enumeration of ® 0 » * 2 * etc * 

the three Green's functions also contain an Iterating enumerat- 
ion that Is signified by a superscript while being separately 
designated by a subscript. Thus, the first Green's function has 
an Iterating enumeration of G°, G^, G^, etc. and similarly for 
the secono Green's function. The third Green's function Iter- 
ates similarly except It begins with the 1 level Iteration. 

The o level iteration In all of the quatlties are referred 
to as the static values since at this Iteration, the frequency 
always equals zero. Additionally, G° and G^ are static 
functions while at further iterations, they become frequency 
dependent functions. The propagation constant 3 has an 
iterating enumeration of 0^, 02 * ^ 3 * etc * There is no 0 q 
since there is no wave propagation on the strlpllne at at " 0 . 


The characteristic impedance Z has an iterating enumeration 

c 

of Z , Z , Z , etc. The frequency and frequency de- 
c o C 1 c 2 

pendent variable (wave number) 



also have an iterating enumeration. The frequency iterates in 

the manner u>^, etc * iterates in the 

fashion k , k . k „ , etc. Since the stripline is speci- 
o,l o ,2 0,3 

fled using cartesian coordinates, the subscripts x, y, and z 
are added to signify particular components of vector quanti- 
ties. As an example, signifies the y-component of the 




17 


vector potential X in the 1 level iteration. Many of the 
quantities are expressed by using Fourier series representa- 
tlon. Each particular Fourier component is assigned a sub- 
script n to mean the nth component. Examples of this nota- 
tion are J and gJ . J denotes the nth Fourier 

ozn 2 n ozn 

component of the static z directed current density while 

denotes the nth Fourier component of the second Green's 

function in the 1 level iteration. Any Fourier components of 

3 A. 

derivatives will be designated as, for example, Two 

frequency dependent parameters which are defined later in the 

theoretical development also sequentially Iterated. They are 

displayed and iterate in the fashion w, , u>, , u . etc. 

ln,i* In, 2 In , 3 

and u>„ , ,0). , etc. Certain values which are deter- 

2n,i 2n,2 2n,3 

mined are not exact values but are rather relative values which 
are subsequently correctly scaled. They are designated using 
a subscript "r" as, for example, A . Where current and 

OZ ft 

charge densities are used to calculate the total currents and 
charges, respectively, the subscript "T" is used to denote 
a total quantity. Thus, 1 ^ means the total 1 level 
current. The del operator is also modified in certain sit- 
uations to V to mean the operation is carried out with 
respect to the transverse components only. In this thesis, 
it is always the x and y components which are transverse. 
The basic field equations initially are developed in the gen- 
eral case with no iterative subscripts. The iterating 
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enumeration subsequently Is applied to these equations while 
developing the proper Green's functions and wave solutions to 
use In the Perturbation-Iteration Method. 

A modification to the calculations Is employed which 
facilitates an easier and faster analysis. With reference to 
Fig. 2.1, all dimensions of the mlcrostrlp structure are 
measured relative to one-half of the width W of the strip line. 
Therefore, the mlcrostrlp is considered 2 units wide, the 

A 

side ground walls extend from -a to a where a « ^ and 

H 

the height of the dielectric substrate equals a where a ■ ^ . 
For certain calculations, the stripline is then located at a 
height y ■ a+ . This modification permits this analysis to 
apply to any size structure with proper scaling applied to the 
dimensions of it. 

2. 2 The Field Equations 

The analysis begins with the formulation of the field 
equations. Since the fields are assumed to be propagating 
down the stripline in the +z direction, they may be represented 
as 


-J0» 

A - A(x, y)e 

± * 

E ■ E(x, y)e 

p » p(x, 

-j8z 

y)e 

-J3z 

* _ -JP* 


-j0z 

$ - $(x, y)e 

H - H(x, y)e 

J « J(x, 

y)e 


The field equations are developed by employing Maxwell’s 


equations 
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7 x t - - I 7 - -J oiy l! 

dt o 

7xS-|-+J - Joie «(y)2 + J 
at o 

7 . $ - U 7-ft - 0 

o 

7 * 5 - e 7* (x(y) I) - p 

o 

We let 

7xH- — 7x7xA--p (77*£-7 2 t)- j«E Q x(y)E 4-J (2.1) 

w 0 U o ° 

Using the Lorentz condition: 

7 *it • -jwu £ x(y) t 
J 00 J 

77 • A ■ -Jwy e 7(x(y)$)« -juiy e [ x(y)74> + $7x(y) ] 

00 00 

and combining with (2.1) results In 

— Ja>e Q £ k ( y > 7«P + #7x(y)] - ~7 2 i£ “ jaie o <(y)f + 5 (2.2) 

o 

Combining (2.2) with S =» -juiit - 7<D gives 

-ju)£ x(y)7<jj -jue *7x(y) — — ? 2 i£ » J + w 2 e o x(y)£-ju»e o K(y)7$ 
0 0 U o 

which can .be written as 

7 2 it ■* v 2 l + 7 2 it = 7 2 J -3 2 A— \i 3 r — tx) 2 p £ K(y)A-juJH e -SVicCy) 
tzt OOO O O 
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2 2 ** 
Letting k Q ■ w M 0 £ q gives the final equation for A : 


7 2 A + («(y) k 2 -S 2 ) X * -u q 3 -juu o e o ^7<(y) 


(2.3a) 


The y component of this equation is 


7 2 A y + (Kr(y)k^-0 2 )A y « -jwp o e o (l->c)4>(a)6(y-a) 


(2.3b) 


Mote that 7*<y) - a y (l~ic) S(y-a) since <(y) undergoes a 

step change from tc to 1 as y crosses the point y “ a. 

3 2 A 

In order that the second derivative term — yields a delta 

3y 3A 

function term, A must be continuous at y * a and ■ ^ 
y ®y 

must have a discontinuity of strength jwvi o e o (l-K) $ at y ■ ct. 
A similar equation for 4 is developed. We begin with 


7* ft 


7 * (tc(y)E) - — 

O 


Substituting in E “ -juiA - 74 yields 


tc(y)7* (-jdjJ-74) » ?‘7ic(y) 


By using the Lorentz condition we obtain 







V 2 * - (tc(y)k 2 - B 2 )* 


E v (l-K)S(y-a) 

•c(y) 


Since the stripline is assumed to he at y ■ a + , the source 

E (1-k) 6 (y-a) 

term K f y) - e - f-> The term ■ * “ ( - y) equals zero 

o o 

at all points y j* a since 5(y - a) « 0 at these points. 
This term may be eliminated from the partial differential 
equation through use of an appropriate boundary condition on 
3<S/3y at y - a, as we demomtrate below. We can replace 
j£’Vic(y) by 


E - -( juiA + li ) 

y 3y y 3y 3y 


and rewrite the equation for 9 in the form 

. .T 3 2 4> . 3 2 $ . 3 2 «1 . 3* 3*(v) . 2, ,,.2 ^ 
+ 1^ + 1F| + 9?* + k (y)lt ° 4 

3 , s 3t / 4 f 3 2 t , 3 2 * L 2, >.2 , 

* ^ <(y) “ai + <(y) ~ + — +< (y)k c < 

3y 3y L 3x 2 3y 2 J & 


- -2- + juiA (k-1) 6 (y-a) 
o y 


since A y is continuous at y ■ a. If we integrate this 
equation once with respect to y about a vanishing small 
interval centered on y-a and take note of the fact that 
p Is assumed to be just above the interface we obtain 
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a 


+ 




a 


jw(K-l) (a) 


This is the required boundary condition on ~ at the 
interface. In each region y <a and y>a we note that 
Vic(y)=0. Hence we can solve for 4 in each region and proper- 
ly join the two solutions using the above boundary conditions. 
Thus the equation for 4 that is used is 


V 2 4 + (K(y)kJ - B 2 )4 . - -£- (2.4a) 

L w 0 

o 

along with the boundary conditions that 4 is continuous 
across the air-dielectric interface and 
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3y 



jw(K-l) Ay (a) 


(2.4b) 


In the lowest order solution 
made continuous at y ■ a. 


is zero so 



is 
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2 . 3 Zero Order Solution 

The Perturbation-Iteration Method commences with the 
evaluation of the static field entitles and Green's functions. 
In the static analysis, u • 0 and therefore S ■ 0. Thus, 
the wave equations reduce to Poisson's equations and are 

7? t “ -liJ(K)i(ya + ) (2.5) 

to o o 


7 t *o 


- 5 (y- a + ) 


( 2 . 6 ) 


At the static level there Is only a D.C. charge and current on 
the microstrip. This translates into the following cond* on: 


3 



a 


z 


Therefore (2.5) reduces to a scalar Poisson equation 


7 t - -Voz<*> s (2 ' 7) 

The two static Green's function, G° and G°, are now 
developed to facilitate the solutions of the equations (2.7) 
and (2.6), respectively. These functions will be developed 
using Fourier series analysis. Only the odd harmonics are 
required since the solution must be even about x = 0. 
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The Green's function g“ is found first and is used to 

solve for J in the integral equation: 
oz 


if G° (x t x',y,y ' 


) J (x',y ') dx' - -2S - 1 
7 OZ ^ \1 

o 


(2.8a) 


obtained by setting equal to p q on the strip. This 

value of A is shown later to be a relative value that can 
02 

be correctly scaled to yield the true value of J . 

OZ 

Since relative values are used, the above Is rewritten as 


21 G° J dx 
I 1 oz,r 

•Jn 


OZ. 


- 1 


(2.8b) 


G“ is chosen to satisfy the equation 


V* G° — <S(x-x , )5(y-y’> 


(2.9) 


with boundary conditions 




are continuous at 


y 


o 


G° is continuous at y ■ y* 

- 0 at y - 0, a. 

G° » 0 at x =» a, -a 
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Assume 



that the solution Is: 

CO 


£ 

1.3 

#2 

£ 


a sinhu ycosu x 
n n n 

-\y 

b e cosu x 
n n 



y i y* 
y*y’ 


In solving for G^, the following Fourier series representation, 
which is employed throughout the entire analysis, is used: 


<X> 


£ 

L-l. 


A cosu x 
n n 


5(x-x* ) 


A 

n 


cosu x* 
n 


a 


Equation (2.9) is integrated across the discontinuity y', 
taking Fourier components only, to give 




cosu x' 
n 


a 


Finally, by satisfying the boundary condition of G° being 
continuous at y' , employing the above, and taking only the 
Fourier components, the Fourier coefficients are solved for: 

a sin u y* * b e 
n n n 
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“‘V* 

-tu b e - tu a coshtu y' 
a n n n n 


cost*) x* 
n 

a 


Thus , 


n 


-a) y* 

n i 

e costu x' 
n 

tu a 
a 


sinhto y* costu x' 
tr n 

tu a 

n 


Since the source points of A are on the stripline, 

02 9 X 

the conversion y' • a will be made. Hence G° is given 

by “ w n a 

e co8tu fl x' 

S sinhtu v costu x y S u 

tu a a n 

n-1,3.. n 

G° ■ , f 

1 sinhtu a costu_x’ -tu y 

n n n . 

e costu x y is 

tu a n ' 


2 

n-1 , 3 . . 


(2.10) 


tu 


nr 


h 2a 


The Green's function is solved for in order to 

solve the equation 


2 J 6°(x,x',y, 


y')P (x* ,y' ) dx' - e * - 1 

'o' o o 


( 2 . 11 ) 


where $ is set equal to ~ on the microstrip. In the 

o 

same manner as A , this value of <J> is nominal and is 
oz o 

used later to properly scale A and J 

r r ' oz,r oz, r. 
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G® must satisfy the equation 


G° -SCx-x’Wy-y’) 


with boundary conditions 


. 3G 2 

g 2 » K(y) If 


are continuous at y * a 


G° is continuous at y' 


G^ “ 0 at y • 0, 00 


G 2 ■ 0 at x * a, -a 


Assume that the solution is: 


E a sinhw y coso> x o S y £ a 
nun 

n-1,3. . 


(b e' 11 * 117 + c e^ 1 * 7 ) cosoj x o.<y?y' 
n n n 


n»l f 3. . 


E d n 

n«l,3. . 


e -Uny cosu) x y's y 

n 


a> 


mr 


n 2a 


As shown before, 
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COSO) x' 
n 


a 


The boundary conditions of continuity of 


sg: 


G 2 and K(y) af" 


at y ■ a, continuity of G° at y ■ y', using the above, 
and taking only Fourier components yields the following equa- 
tions : 


-u a ui at 

a sinhu> a ■ b e n + ce n 
n n n n 


-aj^a ai n a 

koj a coshu a ■ -a) b e +ta c e 
n n n n n n n 


-to y' u y* -w y' 

. _ n , _ ^ n . n 
be + c e - d e 

n n n 


-ta y' -a) y* oj y* costa x' 

cade n + iii b e Q - u>ce n » - 3 — 

n n n n n n 


Once again, by assuming that the source points are on the 
microstrip and setting y' m a, the final solution becomes 


G 


o 

2 


z 

fl® 1 f 3 * * 

CD 

£ 

n“l,3. . 


costa x'sinhcu y 

n rr 

m a(sinhu) a + tccoshw a) 
n n n 

iii a -u y 

cosu x'sinhu ae e 

n n 

oj a(sinhco a + iccoshta a) 
n n n 


COSO) x 

n 


costa x 
n 


0 $ y £ a 

( 2 . 12 ) 


a S y 
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<*» 

n 


Off 

2a 


The numerical solutions for p and J are next performed. 

To accomplish this task, <2. 8) and <2.11) are utilized 

with the respective Green's functions G° and G°. Since 

the potentials are being fixed and measured on the microstrip, 

the observation points of G° and G° are at y ■ a. As 

previously indicated, A Is fixed to the constant value 

oz i r 

of ~ and $ is fixed to the constant value of £ . These 
U Q ° o 

values are nominal and permit easier calculations of J 

oz,r 

and p Q . They basically represent the level of excitation on 
the microstrip and can obviously be varied to any values. 

The values of A and J are properly scaled at 

OZ f IT OZ y IT 

a later point. As Is also shown, the amplitudes of 


both potentials are independent of 8 and Z . The first 

c 


source to be determined 
with (2.8) 


is J . The process begins 
oz,r 


2 f 1 G? J 
I 1 oz,: 


dx' .-2S*£ 


- 1 


« co8u x' sinhw a e 

• y — s — 

I Z— u a 

J 0 n-1,3.. n 


-a) a 

n 


cosoj x J dx' (2.13) 
n oz,r 


o 

Is simplified and approximated to permit an easier evalu- 
ation. We rewrite the expression for G^ as follows : 
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-co a 

-X cosai x'slnhu a e 

z — 2 — 2 — 


n**l |3> * 


to a 

n 


cosai x 
a 


E 

n*l , 3 . . 


-2 to a 

costo x* (1-e n ) 


co a 

n 


COSIO X 

a 


n«l,3. 


. -2 to a 

- [cosw (x-x* ) + costo (xfx' ) ] U-e n ] 
lim xx n 


By applying the identity [ ] 


f e jnu 1 1 + e Ju 

Re z] "V- " 4 Re In — e ' - - 

^ n 2 i _ ju 


n«l y 3« » 


|xn 


tan j 


to G°, the result Is 


G ° ■ ~ 4^ ln i fan l tan ^■t 2efx ' I] 


4a 


4a 1 


- ln< 
8ir 


(*• , ira , 

cosh + cos — 

a 


-A 


. ira 

cosh ~ — cos 


9 ^(X-X')ircosh “ + COS 
s -“(x-x' jJ^cosh ~ - cos ~(xfx' )_ 


G? + G° 
la lb 


The series Identity 


2 ,4 

lnjtanu | « ln|x| + + . 


is applied to G^ a » the dominant terms are retained, and the 


< ■ V, ■* ivnOim ; *■' . ' 1 
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result is 


G 0 - - In ~ 

la 2 tt 4a 


fa In | (x 2 -x ,Z ) 


“ 9^7 < x2+x ' 2 > 


iso? <ir)*<»Wx' 2 4«' 4 ) 


The series identities 


2 4 

X X 

costuc ■ 1 + Yi" + + • ■ • 


cosx - 


2 4 

i x . X . 

1 “ 2! + 4! + * * * 


are applied to G^, the dominant terms are retained, and the 
result is 


G °b" ~k ln ln (t<2o) 2 + (*-x’) 2 ][(2a> 2 + (***') 2 ]j 


+ 757 <77 )4 K*-*’> 4 + C=«-*’> 4 ] 


r 2 /TTv2 2 , tt v4,. .2,. 2, ( 2 S 

+ 1 U W ' 77 ( u> <2a) > (x **' > 


2 /« »2 1 »4 

' 37 W (2 “> + 367 <4l> (2o) 


On the stripline, and ^ are very small provided a is 

chosen to be greater than or equal to 10. The two functions 
o o 

G^ a and G^ are recombined and, after employing the fact 


ii 
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that terms containing (^) 4 are negligibly small and thue dis- 
carded, the final result is: 



f( 2 a) 2 + (x-x’) 2 ] 


[ ( 2 a) 2 + (x + a') 2 ] } + 


1 ( J L) 2 
3 Ua ; 



(2.14) 


The current density J qz r is approximated by a series 
of the form 


, (,.) . V V' 2 * V s v ’ 6 

yfrr^ 


This form gives J „ „ the required singular behavior at 

oz,r * 

x " + 1. The coefficients 1,1,, I., and I, are solved 
■* ox/ j 

for to determine the current density J . By substituting 

OZ | V 

and J OZ(r back into (2.13), the result is 



- -In (( 2 a ) 2 + (x-x' 


) 2 ] [ ( 2 a ) 2 + (x +x’ 




Point matching on the microstrip is used to solve for the 

four unknown current density coefficients by performing the 

113 

integrations, sequentially letting x ■ ^ ^ * " 4 * atu * 
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and chan using Che four equations Co solve for Che four unknown 
coefficients. The individual integrals are solved utilizing the 
implicit transformations: 

x • sin0 
x' • sine' 

x 2 -x' 2 ■ y (cos20' - cos20) 
dx' - c.os0'd0* 


lnt-jCcoa©' - cos0)] 


20080000800 ' 


ZcosnB 

n^2.. n 


- In4 


Thus we find that 


- II iajT i ° <* 2 - *' 2 > to’ - la ln2 

iiirv 

■ ^ Jo 


yjV- x' 2 


2 ln( x 2 ~' X ill dx' - ~ (2x 2 + ln4 - 1) 

>/r^”x* 2 


W - x' l ) dx> + ln4 , 1 

4tt In I 8 4 8 


2I 3 f x* 

'Wo 




6 ln(x 2 - x' 2 ) 


x' 2 
I 

64 v 3 


dx' 


r£( ^ x 6 + 4x 4 + 6x 2 - + 51n4) 
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^> 2c2 “) 2l of l -7= 

J o V 1 '*" 

1 2 /TT.2 r - .2 t f 1 x' 2 
2i I W (2a) h r=-= 

J Jl ~ 

0 > 

l 2 .IT ,2,„ ,2 r r 1 x' 4 
2? 3 ( «S> (2a) *2 / ^== 

{ 


i ! $;>*»')% f l -j===. 

{ i 1 '*' 


- i <*> 2 c 2 “>\ 


dx* - 


dx' - 


i <£> 2 < 2 »> 2i i 


I 4> 2 < 2 “> 2i 2 


■**' - ! 'ife )2(2al2l 3 


The integral terms 


-faj In I (2a) 2 + (x - x’) 2 ] [(2a) 2 + (x + x’) 2 ] 
fl + I,Jt' 2 + I,x’ 4 + I,x ,6 1 , 

L >/r^ J 


cannot be solved in a closed form fashion and, thus are numer- 
ically evaluated using a Simpson’s Rule approximation. The 
substitution x’ ■ sin0’ is employed to avoid a singularity in 
the calculation at x* - 1. 

Once the current density coefficients have been solved 
for, the total current on the microstrip I qT r is evaluated 
by using the integral on the following page: 


\ 
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and Che substitution x' - sine. This yields the result: 


W" tI o + l I l + l I 2 + 




The charge density p o and subsequent total charge q QT are 

solved for in a similar manner as J _ and I _ . The evalua- 

oz t r oi i r 

tion begins with (2.11) 


1 

G? p dx f - e * - 1 (2.15) 

L O O O 

G^ is simplified by expressing it as a function of G°. 



. cosu x'sinhw acosw x 
o _ X™* n n n 

2 




at a(sinhu a + tccoshw a) 

n»l,3. .1 n n 


i- k ♦ ti 

L n“T,3.. 

/ K + 1 

\sinhoj a + KCoshw a 
' n n 


-1- cosui x'sinhw a 


- 2 e 


-w as I 

i ICOSU) X 

) n J 
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The charge density p o G£") is approximated by the series 


»„<*’> - 


"o + V 


- 2 + 




. 4 r 6 

4 2 * + 4 

“72 


which is similar to that used for the current density. Thus, 
in the same manner as before, the four unknown charge density 
coefficients are solved for by substituting the expressions 
for p and into (2.15), performing the integrations, and 

O £ 

solving the four equations obtained by point matching at x ■ 
113 

j ~ , and 1 for the unknown charge density coefficients 

q o 1 q l ’ q 2 ’ and q 3* 1116 inte Srals 



dx» 


are performed in the same way as before using the approxima- 
tion of G° and the Simpson's Rule approximation. The remaining 
integrals 


1 

f y 

_4_ 

sinhw ot 
n 

r < +• 1 -wo 

-2e 

0 n=l,3. 

nn 

• 

K + 1 

sinhtjj cx + k coshui a « 

n n 


COSU) X 

n 


.2 


.*4 -L. .. „« 6 


COSO) X 

n 


, q Q + q l X + q 2 X + q 3 X 


dx' 


y/l - x' 


are carried out using a Simpson's Rule approximation and the 
substitution x' ■ sin0' to avoid the singularity at x' * 1. 
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The total charge on the microstrip Is determined by 
using the Integral: 


q oT 



1 

P 0 (x’) dx' 


0 


and the substitution x’ - sind*. This yields 


^ « 5 

q oT ‘ T(, Io + 2 q l + I q 2 + 16 q 3 ) 


With the results of I qT r and q QT , the static Z cq ^ 

k CJ . is detetmined. At id ■ 0 , 0 0. Therefore, 0. is in- 

eff,o o a 

corporated into the analysis with the frequency being set equal 
to zero at the conclusion of the static case. Beginning with 
the equation 


E * -jw.X - V$ 
o J 1 o o 

E - “j wA + j0 $ 

02 OZ 10 

On the strip, 

E oz ■ 0 ■ * J " A oz + 1 * 1*0 ’ ‘ i “ L o I oT + 1 9 1 "c 1 C2,16) 

o 

In this equation, L = A /I „ is the static inductance per 

0 OZ 01 
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unit length of line and C q ■ q o ^/® 0 is the static capacitance 
per unit length of line. Using C2.8) 



it is evident that A and J are directly proportional to 

o z oz 

A and J , respectively. Thus, the value of L can also 
oz,r oz,r r o 

be determined from L » A / I _ . This illustrates why the 

o oz»r oi, r 

values of 8. and Z , as yet to be determined, are indepen*'- 

X CO 

dent of the nominal value of A . The value of C is also 

oz,r o 

independent of the level of * by the same reasoning. From the 
continuity equation 


* 3 


3J 


ox 


3x 


- 


1 oz 


-j Vo 


Integrating over the width of the microstrip yields 


VoT " “l q oT <2 - 17) 

since J *0 across the strip. Combining (2.16) and (2.17) 

OX 

finally results in 

8, - uX C~ (.2.18) 

1 11 O O 


This is the value of the propagation constant in the zero 
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frequency limit and is also used as the approximate value in 
the 1 level iteration. As can he seen, 0 at Wj- 0 which is 
consistent with all of the previous theory. 

A value of 0 is determined at o>j- 0 by defining it 

as follows: 



(2.19) 


The final parameter to solve for is the characteristic 

impedance Z of the stripline. This is easily determined by 
co 

the following: 


Z 

co 



VoT 

“l C o 


oT 



( 2 . 20 ) 


The current and charge density coefficients are now 
scaled correctly as they are needed, properly scaled, for fu- 
ture calculations. Equation (2.16) 


u)A 


oz 


e i»0 


properly relates the levels of the potentials. Let 


A - KA 
oz oz,r 
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where K is the constant of proportionality which correctly 

scales A . As stated previously, 
oz 

A oz,r " p o 

* - - 7 - 

o e 

o 


Thus 


K 



WU G 

o o 


which, from (2.8), is also the constant of proportionality be- 
tween J and J . Therefore, the current density and, fur- 
oz oz ,r 

thermore, total current I qT is properly scaled by the relations 


J 

oz 


\l L C 
y o o 


j 

oz,r 




oT,r 


(2.21a) 

(2.21b) 


2.4 Higher Order Solutions 

The analysis continues with the calculation of $2 anc * Z^ 
at a finite frequency <o^. The value of should be low enough 
to use the static solution as a good first approximation be- 
fore calculating a better value of 8. The procedure continues 
with the calculation of the Fourier coefficients of p and 
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J . This Is done using 
oz 




(x')cosu x'dx' 
n 



(x'Jcosw^x’dx 1 


0 


(2.22a) 

(2.22b) 


along with the substitution x' - sine*. These integrals are 

solved numerically using a Simpson's Rule approximation. A 

sufficient amount of coefficients should be calculated until 

p and J are very nearly equal to zero, 
on on 

The Fourier representation of A qz and 4> o are solved for 
at any point in the cross-section of the stripline. A qz 
satisfies (2.8) 


7 t A z 

t oz 


“ lJ o J oz <x),S(y " <x) 


with the boundary conditions 


9A 

02 

A , are continuous at y ■ a • 

oz 3y 


A * 0 at y - 0,® . 

oz 


A * 0 at x » a, -a . 
oz ’ 


A Fourier representation of A is used. Assume 

oz 
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L A sinhw y cosw x 
n n n 

n»l , 3 • • 

oz „ -« y 

E B e cosu x 
n n 

n- 1 9 3 • • 


y ( a 


y S a 


The coefficients, A and B , are solved for by invoking 

n n 

continuity of A at y *a and integrating (2.8) over a small 
oz 

region centered on y ■ a , thus 


A sinhw a * Re 
an n 


-ui at 
n 


3A 


ozn 


3y 


-u a 

-oj B e n - n A coshw a * -U -J 
nn n n n o on 


A " — 
n w_ 


y —W a 

o n . 
e J 


on 


n w 


sinhw ot J 
n on 


Therefore, the Fourier representation of A qz is 


oz 


V- u o “V 

<L ~ e J on alnh “n° C ° 8 V 

n»l,3. . n 

A Ho _ <V 

/ — sinhw a J e cosw x 

Lm* w _ n on n 

n- 1 , 3.. n 


y $ a 


y 5e a 


(2.23) 
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$ is found using (2.11) 


„2 * 
7. 9 
t o 


- 7 1 fi(y-y') 

o 


with the boundary conditions 


3« 


9 • <(y) t — continuous at y ■ a 

o 


9 q continuous at y ■ y' 


4> » 0 at y • 0, ® 

o 


"0 at y ■ a, -a 


Assume 


Z a sinhat y cusw x 0 SV{ a 

n n n 

n-1 ( 3*. 


9 “ 


°° ^ -oi y w y 

(be n + ce n ) costu x a s y s y' 
J-f n n n 

n»l, 3. « 


.y 


E v 2 

n«l,3.. 


COSU) X 

n 


y's y 


The equations used to solve for the unknown coefficients emanate 


from the continuity of <t at y = ct, the continuity of 
3$ 


< (y) ~~L at y » at, the continuity of 9 q at y = y* , and the 
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Integration of (2.11) over y. Thus, 


-<u a (j a 

a sinhw a ■ b e n +ce n 
n n n n 


-a> a w n a 

tea cit co BO) a ■ - b in e + c w e 
n n n n n n n 


y' <*> y' • -w y' 

be + c e ■ d e 

n n n 


-w y -u y' w y’ 

-dwe . b u e +cwe ) 

n n v n n n n 



After solving the equations for the coefficients and letting the 
source point y' « a , the result becomes 


p slnhu y 
on n' 


„ . - E at (sinheo a + iccoshai a) 
u"i y J • • o n Q n 

*o « p e tJn ° l sinh<i) a e n 

on n 

e to (slnhu) a + Kcoshw a) 
n«l,J.. on n n 


COSUl X 

n 


COSU) X 

n 


y s a 


y % a 


(2.24) 


The analysis continues with the next iteration of field 
entities. The frequency is raised to and the propagation 

constant is approximated using the static solution of 3^. The 
Green’s functions are also redeveloped as they are now chosen to 
he frequency dependent. All of the following theoretical develop- 
ment will be enumerated as the 1 level iteration. It should be 
noted that all further iteration processes will be carried out in 
the exact same format. 


i.YA'Iwf? , W&Sd£j2 
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In che zero frequency limit A and A are both zero. At 

ox oy 

this Is no longer true, but nevertheless the first correc- 
tions for A^ and A^ are small. It Is possible to find both 
A^ ^ and A^ approximately by using the solutions for A qz and 

4 . We can then find better approximations to A and 4, which 
o z 

are labelled A , and 4,. We can then Iterate again to obtain 
lz 1 

Improved solutions for A^ and Aj y and then for and 

After a few Iterations the solutions converge to stable values, at 
which point the frequency u Is increased to The solutions 

obtained at are used In a similar Iteration process to find 

the converged solution at w ■ 

A^y Is the first field entity to be found. From (2.3b), A^ 
must satisfy 


I** - (*i - v 


Join e 4 
o o o 3y 


(2.25) 


where 


. . (K-DSCy-a) 


Note that the known potential 4 q is used on the right hand side. 
The solution for A^ y is 




[-jwy e (k- 1) 4 (a)] dx* 
oo o 
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The Green's function will be chosen to satisfy the equation 

G* - - 6(x-x') 6 (y-y’) 
with the boundary conditions 


i 30 $ 

G-, - — are continuous at y ■ 
3 3y 


3 Gt 

■ 0 at y ■ 0 

3y 


G^ "0 at x • a, -a 


(2.26) 


Assume 


°3- 


OD 

y 

A 

cost!) x' 
n 

coehw 2n,l : 

n-1,3 

A 

n 

• • 

a 

00 

z 

n-1,3 

B n 
• • 

cosu X* 
n 

■w, 

ln,l 

■ e 

a 


COSU X 

n 


y S a 


y i a 


where 


- <k^ .+ S, <u. \ 

2n,l N n o,l 1 ln,X > 


2 2 2 

</ - k z . + e: 

n o,l 1 


The unknown coefficients are determined by using equations satis- 
fying continuity of at y » a and the integral of (2.26) over 

a small interval along y centered on y “ a . Thus 
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A coshu) 
n 


2n,l 


B e 
a 


-o>, .a 
ln,l 


“w, -a 
n In • X 
- 0), .Be ’ 
ln,l n 


“2n,l A n alnh V.,l a ' ' 1 


Hence, 


A - 

n 


“2n.l alnh "2n,l ct + “ln,l coah “2n,l“ 


B - 


C °a htlJ 2n.l ae 


<D. ,a 
In, I 


o)„ , sinhoi. -a + w, ,cosho)„ .a 

2n,l 2n,l ln,l 2n,l 


and the solution is 


L - 

n-1,3. . a (o> 


cosw x'coshw. ,y 
n 2n,l 


E T 

n-1,3.. a W 


„ ..sinhu),. ,a + w, , coshoi- ,a) 
ztif x zn f x in fX 6ii|i 

ai. .a -w. . y 

cosoi x'coshoi. .a e ’ e n ’ 
n zn.i 

2n,l 8lnhu 2n,l“ + “ln,l COa "2n,l a) 


cosoi x y 
n 


cosoi x 
n 


where 


0) 


ln,l 



0) 


2n,l 




+ 6 ' 


S a 

( 2 . 27 ) 
i a 
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The solution for A, at y • a is 

ly 


* ly “2n.l slnh “2n.l a + "ln.l COBh “2n,l“ 


COSO) X 

n 


The quantity 


3A^ 

3y 


can now be determined and is 


y - a 


3Aj 

3y 


+ ■ £ * “ln,l A l)4 o,c ° 0 V 

y • a n-1,3.. 


3 Vo e o M ln.l (,: ~ 1>cosh "2n.l 3 *ou^> 


l^o o ln.l 2n,l on' 

. . “2n.l slnh “2n,l a + “ln.l C “ h, Vi 


(2.28) 


coso) x (2.29) 
a n 


The quantity $^(x,a ) is determined in terms of A lx . On the 
microstrip, E^ x ■ 0. Thus, 


3« 

E u ■ - W - 0 


Integrating over x yields 


$ 1 (x,a + ) » -jo A^dx* + ^(0,0) 


0 x x* 


- / 

0 0 


3A. 

^r# dx” dx' + ^(O.a) 


(2.30) 


The term 4^(0, a) represents the D.C. value of which is maintained 

at — , as used in the static solution, for this and all other 
o 


..n.i/.-i-.t, - 
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7 • it 


3A 
2 

3x 



AA 1* 


function 


found using the 

3A 

3A 

3A„ 3A 

ji 

z 

£ + — Z _ 

3y 

32 

3x 3y 


-Jo»V o e o ic(y)« 


At y ■ a , the solution for 


3x 


becomes 


- J8 l A o* - w* ' 


3A, 


(2.31) 


After representing (2.31) as a Fourier series, substituting it into 
(2.30), and performing the integrations, 4^ is expressed as 


. _____ J(tj« wn. 

> ’2 2 a6 X A o.n - W** * 3 "lVo*on> 

n-1,3.. 


3A, 


(cosw x - 1) + 
n 


(2.32) 


The next step is to solve for the new charge density p^. An 
equation must first be developed to accomodate this task. 
Consider a Green's function G^ satisfying 


I7 t " (3 i " K(y)k o i )] G 2 " 


(2.33) 


with boundary conditions 


1 9G 2 

, tc(y) are continuous at y 
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is continuous at y ■ y' 
Gj - 0 at y “ 0, ® 

G^ ■ 0 at x ■ a, -a 


Assume 


cosu x 

/ \ A. — — sinhoi 


n»l,3. . 


n a 


2n,l y 


cosw x 
n 


0 $ y $ a 


J ln,l J 


oj. ,y cosu» x 
Xu f x \ n 


£ (3 n e iU ’* + C_e *“•* ) r — coaw n x 


a 


a g y s y' 


n™l , 3 • • 


“ cos'*; x 1 -to. .y 
V* _ n ln.l 

) D — e cosw_x 


n“l - 1 . . 


n a 


n 


y’ < y 


where 


a; 


ln.l * 




to 


2n,l 




. + 

o.l 


The equations employed to solve for the unknown coefficients 
use the boundary condition that.G. is continuous at y ■ a, the 

3g£ 

boundary condition that tc(y) t — is continuous at y - ct, the 

dy 

boundary condition that is continuous at y " y* , and the integral 
of (2.33) over y. Thus, the equations are 
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-w. ,a to. ^ 

A sinhw, .a - B e An,A + C e An,A 
n Zn , x n n 




-W- -a cu, .a 

- u .B e, 111,1 + «. e ln ' 1 
ln,l n ln.,1 a 


,y' to. ,y* -to, ,y' 

B e 111,1 + C e 1x1,1 - D e 1x1,1 

n an 


- w- -D e 
ln,l n 


-to- ,y' 
In, l 7 


-co, .y' til, ,y* 

“m,lV ln • 1 + • ) 

* - 1 


After solving for the coefficients and letting the source point 
y' ■ a, the solution for becomes 


cosco^x* sinho)2 n 


n-ft.. a(u ln,l 8lnhw 2n,l 0l+Kw 2n,i 


„ . — r coso) x o 5 y 

coshco„_ , a) n 7 


2n,l 


w. ,a -co. 


, . . 'ln.r “ln.l 7 

cosco x sinhco- , ae e 

Zft|i 


(2.34) 


n-ft.. aCu ln,l Sinha, 2n.l ct + K u 2n,l COshaJ 2n,l 


Next consider the following: 


[^2 * <“n ' K( * )k o,l + 6 J>] =2n - 


—7 cosco x a $ y 
a) n 7 


In 


S 2n[^2- <“a - ^ k l,l + 6 1>] *2n 
dy j^ln TT 2 - G 2n + =2n T* 
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Integrating over y yields 


8g L 

. 

— — * 

In 3y 

u 2u 3y 

, 8G 2n 

r l 3 * 

In 3y 

" ^2n 3y 


In 


In 


’+ - - »m<‘ + > + G k<“ + > t 2 

a o 

a 

- 0 

0 


(2.35) 


(2.36) 


Multiplying (.2.36) by ic and adding to (2.35) gives 


1 . . 3t ln 

G 2n ST - 


- * - + G 2n ( “ +) -t 

a o 


Using (2.4) 


C 2n<“ + > Tf - V“ + > - >1 (1 -> G L A lyn ( ° +) 


(2.37) 


Summing the Fourier components (2.37) results in (y ■ a ) 


t ' ° 

n»l, 3. . 


_1 M ln 

°2n T- - 

0 n-1,3.. 


«> 

E [*ln - >l <1 -’ t)G 2u A lyn] c03 V 


(2.38) 


Equation (2.38) is thus equivalent to 


7 / G ^ 


dx’ 


*1~ 


00 1 

J^u-K) £ 0 2n A lyn «=V 

n-1, 3. . 


(2.39) 


Equation (2.39) is used to solve for the charge density p^« 

The procedure in carrying out this assignment begins by substituting 
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(2.32) for is 

nent in (2.28). is 


substituted by the equivalent Fourier compo- 
substituted, letting y - o, by 



slnto2 n-1 » 

“ln,l alnh “2n,l a + KW 2n, l coshu 2n, 1° 


As done previously, is represented by the form 


p i Cx ' ) 


q „ + i ^' 2 + q 2 *' 4 + < l 3 1 ' 6 


Vi - x ' 2 


The integral on the left side of (2.39) is carried out by 
making a substitution for G^. This substitution allows a rapid and 
accurate convergence of the Integral. The substitution is as 
follows; 

2 J G^dx’ -2 J G^Pjdx* + 2 J (G* - G°) P] dx' 

0 0 0 

a 2 / G 5f*i d ** +2 f E (G 2n “ G 2n )p l dx?c0aw n x 

n n n*l,3.. 

0 0 

* 2 J G 2 P + 

0 

f rcos(i> n x f sinha) 2n 

2 in-l1..L a< “ln,l 3lllh “2n,l C ‘ + K “2n,l coah “2a,l°‘ ) 


0 

cosai x'sinhw a 1 

n n | 

suu (sinhcu a + Kcoshai a) I p l' 
n n n •* 


dx'cosai x 
n 


f 1 - f 1 

21 G° Pl dx' +2 ^ C a I pj^cosw^x'dx'cosw^ 

J n-l73.. J 


x (2.4C) 


54 


The Integrals Involved In solving for 2 




dx' are identi- 


cal as the closed form approximation previously solved for in the 

static solution. As n gets large, the value of quickly approaches 

zero which facilitates a rapid convergence of the Fourier series 

summation. The integrals involved in the summation are evaluated 

using a Simpson's Rule approximation. 

The charge density coefficients are determined sing (2.40) in 

113 

conjunction with (2.39), point matching at x ■ ^ , j ^ , and 1, 
and solving the four resulting equations for the coefficients. 

The current density J^ z is determined by first solving for a 
relative current density J, and then properly scaling it to give 

J.2 1 XT 

J, . J, has the familiar form 
lz lz,r 




I + I.x’ 2 * + I-x’ 4 + I.x' 6 
o 1 Z 3 




From (2.3), 


hz.r ■ - "o J lz,: 


where A, is the relative value of A- which is also correctly 
lz,r lz 

scaled at a later point. Therefore, the solution of J lz r satisfies 
the equation 


2 f G*J dx' = 

J 1 l*.r u 


(2.41) 
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The Green's function Is next solved for. must satisfy 
the equation 


17* - (S l - gJ - - fi(x-x')6(y-y'> 


(2.42) 


with boundary conditions 


i 3G i 

G^ and are continuous at y 


G^ is continuous at y • y* 


g| ■ 0 at y ■ 0, 


Gj ■ 0 at x ■ a, -a 


Assume that (y ' - a) 


COSO) x' 
n 


' n a 
i n-1, 3. . 

G l - 

“ cosut x' -uj, .y 

V B a-e ^ 

a a 

n-1,3.. 


sinhai- . y cosu x 
2n,r' n 


COStl) x 

n 


y ( a 


y % a 


r 2 2 2 n 2 2 

1D 1 . - k . + B, 1 = \/ w - + 6, 

ln,l y n o,l 1 2n,l \ n o,l 1 


where 
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The unknown coefficients, A q and B q , are found using equations 
satisfying the houndary condition that is continuous at y ■ a 
and the Integral of (2.42) over y. Thus, 


V lntu “2n.l“ ' V 


~“ln,l a 


- oj- .Be 
ln,l n 




“2n,l A n co8h “2n,l° 


1 


The coefficients are readily found to be 


“ln.l Sinh “ 2 n,l a + u 2n, 1 c08 h“ > 2n,l a 


sinhu>„ , e 
zn^i 




w ln,l stIlhw 2n,l a + w 2n , l cos8w 2n , l a 


Therefore, is expressed completely as follows: 


cosu x'sinhw- ,y 

n 2n,l 7 

a(a)- ,sinhw_ + a, .coshw- .a) cosa) tl x y * a 
n®lj3* » Xn*X 2ii| 1 2n t X 2n*X 


oj 

E 


® cosu x* sinhw„_ ,ae 

V 

a(w. _ 

n°l,3. . 


2n_j_l 


“ ta ’ 1 V “ ln > 1 


(2.43) 


( “ln,l 8lnh “2n.l a + “2n,l cosh "2n > l“ ) 


cosu x y % a 
n J 


In (2.43) we let y * a and substitute this into (2.41). The 


integral 2 


Jn 


is carried out by making a substitution 


1 *'0 

for G^ that gives a simpler evaluation. This substitution is as 




follows : 
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2 J ihz.T^' - 2 f G ? J U,r d *' + 2 j ' G I - °J> J a.r ,ta ’ 

0 °ri 0 « 1 

■ 4 a lhz,^' * 2 t, f < - G L> J l*,r‘ 1 *' 

J n-1,3. . J 


G I J lz.r d *' * 


( 

o 1 

«- / Isinho)- , a 

2 I I 2n .1 

I a (a). .sinhw- -a + oj_ .coahaj- , 
n-1,3. +Jq L ln,l Zn,l 2n,l 2n,l 


sinhoj ae 
n 


-in a 

n 


iw n 


‘ 

J- COS!) X 
X Z y IT Q 


'dx'cosw x 
a 


rl a rl 

2 / G?J, dx* + 2 Y* / D J, cosw x'dx' cosat x 
I 1 lz,r I n lz,r n n 

J n-1 , 3 . . J 


0 

1 


The 


fl 

integrals Involved in solving 2 G I J 1*., c dx’ 

Jo z,r 


are Identical 


as the closed form approximation previously solved for in the static 
solution. As n gets large, the value of quickly approaches zero 
which enhances the rapid convergence of the Fourier series summa- 
tion. The integrals involved in the summation are solved by using a 
Simpson's Rule approximation. 

The value of A^ z r is chosen to equal displayed in 

(2.32). This selection of A. accomodates a simple solution to 

XZ f IT 

J, and is shown to be a completely arbitrary specification. 
iz,r 

The solution of J, is carried out by performing the affor- 

XZ p IT 

mentioned substitutions, creating four equations by point matching 
113 

at x « ^ , j » ^ , and 1, and and finally solving these equations 
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for Che four unknown coefficients. Thus, JL Is now determined. 

LZ | X 

The new and more accurate propagation constant 6 2 and charac- 
teristic Impedance are next determined. The first step in 
accomplishing this is to calculate the total charge and total 

relative current 1 1T on the microstrip. These tasks are carried 

li t r 

out using the integrals 


* 2 j P x (x' )dx' - ir(q Q + \ ^ + g 
0 

- 21 J. (x^dx’ - tr(I +i I, 

,r I lz,r o 2 1 8 


(l 2 + 16 q 3 } 


X 2 + V 


in conjunction with the substitution x' » sin6'. On the 
stripline! 


E lz " + j0 2*l " 0 


( 2 . 44 ) 


As has been previously defined, 


A, - KA, 
lz lz,r 


(2.45) 


where K is the scaling factor. By virtue of (2.41), it is seen that 


J. - KJ. 
lz lz,r 


(2.46) 


Combining (2.44) with (2.45) yields 


Using the previously chosen value 


A 


lz,r 




1 


(2.47) hecomes 



(2.48) 


In the same fashion as (2.17), upon integrating the continuity 
equation over x, it is required that 


0 2 X 1T " W l q lT 


(2.49) 


Substituting (2.46) and (2.48) into (2.49) finally results in 


8 


2 



(2.50) 


This la the new and more accurate value of the propagation constant. 
Using 8 2 , a value of ic eff ^ is solved: 



The final entity calculated is Z £l at the given frequency u^. 
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Z ia defined as the ratio of the line integral of along x ■ 0, 
from 0 tor a along y, divided by the total z directed current I 1T on 
the microstrip, this line integral, which is the total potential 
difference between the ground plane and the microstiip, is found 
from 


IT 


■ J t • dy ■ -J * (- ja)^ - • dy 

r r 

-J (Jw! A ly + jj- > d y - *i<«) + i^J A ly d y 


The total current is solved using the previously determined ^ 


and I 1T>r . 


^2 

I 1T " KI lT,r “ ^ I lT,r 


Thus, 


'cl 


_1T 

C 1T 


(^(^(a) + j 




VlT.r 


(2.51) 


The expression for 


L 


A^ dy is easily evaluated using the equation 


for A^ and letting 0 $ y ( a. 

A property that should be noted is the independence of the 

final results of 311(1 Z C 1 on the assigned value to A^ z This 

independence is evidenced in (2.47) where any assigned value of 

A, would just adjust the value of K. The chosen value of 
iz , r -i- 
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obviously simplifies Che subsequent calculations. Once the value of 
&2 Is determined, the scaling constant K can then be easily calcu- 
lated. 

Before the next iteration process begins for calculations of 
new values of 6 and at a new or same frequency, certain proce- 
dures are taken. The values of the current density coefficients 
are properly scaled by using (2.46) and (2.48). This yields 

®2 

J lz " <!)^ J lz,r 


The Fourier coefficients of and are calculated using the 
equations 


^lzn 

p ln 



ji z (x f Jcosu^x'dx* 

, (x')cosw x'dx' 

1 n 


along with the substitution x' - sinb* and a Simpson's Rule approx- 
imation. An adequate quantity of these coefficients should be 
calculated for proper convergence of future calculations using them. 

An equation for is determined using the properly scaled 
version of (2.41). We begin with 




1 

G?\J dx* 
1 lz 


CO 

V* u G. J, cosu) X 
o In In n 

n-1,3. . 
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Therefore, at y ■ a , 


. . . & Vln 8lnh “2n.l 

A i* ( *>°> ■ E 


w, .sinhw* ,a + a», ^oshw. ,, 

n-1,3. . ln.l 2n,l zn,l 2n,l 


— co8u) x (2.52) 
a n 


where 


(it 


In , 1 





Kk o,l + S 1 


The Fourier series of is next determined. Stating (2.39) In 
complete Fourier series notation. 


»l ( *- o) ' T E G 2n s ln“ 8 V ' E G 2n A lyn c ° 8 V 

0 n"l,3. . n _ l,3. . 

(2.53) 

Upon substituting the previously solved Fourier coefficients into 
(2.53), the final result is 


^(x.a) 


« (p. - jat.e (tc-l)A. ) 

In 1 o lvn 

Z-ii e (at. -sinhat. ,a + -coshaj„ .a) 

n«l,3. . o ln,l 2n,l 2n,l 2n,l 


COSbJ X 

n 


(2.54) 


At this point in the analysis, the 1 level iteration has been 
completed. All entities which are required for the next iteration 
have been solved for. The next iteration may fulfill one of two 
tasks: it may improve the values of 6 and by solving all of the 

equations again with (J- replacing g, and keeping the frequency con- 

Z 4 . 

stant at or it may solve for new values of 0 and Z £ by increasing 
the frequency to a higher and using the approximate &2 at this 
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new frequency* The Iteration repeats at the location of the 
solution of with all enumerated Iteration levels Increased 

by 1 and continues until the Fourier solution of the new $2 
function Is found. This Iteration process may continue Infinitely 
with the frequency continuously Increasing. Figure 2.2 summarizes 
,-».id depicts the order of solution of the field entitles In the 
Perturbation-Iteration Theory. 

This Is the conclusion of the Perturbation-Iteration Method. 
Chapter 3 presents data of various entities calculated from applying 
the theory to different structures. It also compares certain data 
to established theoretical and experimental data to exhibit the 
validity and accui acy of this technique. 























CHAPTER 3 


CALCULATIONS AND RESULTS 

The Perturbation-Iteration Method was used to find the 

effective dielectric constant and characteristic Impedance for 

W 

various combinations of K and dimension ratios — over a range 

of frequencies. In each Instance, the dimension A equals 10W or 

10H, whichever Is larger and H ■ 1mm. This gives validity to the 

assumption of using the side ground walls. The data Is organized on 

the following in cables which list the computed values of K e ff and 

Z c « Each table represents a specified combination of k and 

ratios of the stripline structure over a range of frequencies. 

Graphs of ^ and Z c are displays ’ to interpolate data between 

W 

calculated points and are organized with all combinations of the 
ratios at each value of ic shown on each graph. These graphs are 
labelled as Fig's. 2.3 - 2.14. 

The evaluation of the accuracy of the numerical procedure used 
In this thesis Is difficult because of the lack of good data to 
compare the results with. Xuester and Chang have noted that there 
are large variations in computed values for the effective dielectric 
constant as presented by different authors {18] . This lack of data 
has thus limited the verification of the accuracy of the numerical 
procedure used to that of a few special cases. 

Figure 2.15 illustrates a graph which compares data calculated 
from the Perturbation-Iteration Method with data taken from the 
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following sources: a fullwave analysis by T. Itoh and R. Mittra [14]. 

a dispersion model analysis by W. J. Cetsinger [10], an Improved 
analysis of [10] By T. C. Edwards and R. P. Owens [19J, and experi- 
mental data 119]. The Perturfiation-Tteration Method compares very 

favorably with all of this previous work for the special case of 

W 

k • 10. 15 and three values of g . In the Table accompanying 
the figure 2W and H are given in millimeters. In Fig. 2.11 the 
results are compared with computations made by Jansen for k ■ 9*9. 
Jansen's results were read from the graphical data presented and 
multiplied by 1.01 to make the results equivalent to chose obtained 
for < - 10 using the Perturbation-Iteration Method. The two sets 
of results compare very well for g ■ 1. For g ■ 0.5 and 0.25 
Jansen's results appear to be approximately 1% higher (this con- 
clusion is subject to the error in reading data from the curves 
and may not * e correct). The effective dielectric constant at 

zero frequency, as obtained from Wheeler's approximate theory, is 

w 

also shown for g - 0.5 and 0.25 in Fig. 2.11. These values lie 
below those calculated using the Perturbation-Iteration Method and 
seem to suggest that the latter are most likely very accurate. 

Examination of the tabulated data shows that for wide strips 
and large dielectric constants at higher frequencies the numerical 
procedure used converges to a value for the effective dielectric 
constant that is larger than that of the substrate. This is a 
clear indication that a 4 term expansion of the charge and current 
distribution along with poixxt matching at 4 points on the strip 
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in the interval 0 to 1 is inadequate for wide strips and large 
dielectric constants. The region where the results become inaccu- 
rate are mostly in regions outside of those of practical Interest. 
For example with a substrate having k - 10 a substrate thickness 
of 0.25mm would be appropriate in which case the 5 Ghz data in 
Fig. 2.11 becomes 20 Ghz data (the frequency scale must be multi- 
plied by 4 because of the reduced height). 

It can be seen that some of the dispersion curves for the 
larger values of ^ in Fig’s. 2.5, 2.7, 2.9, 2.11, and 2.14 show 
a second inflection point. This marks the region where the numer- 
ical solution begins to become inaccurate. The boundary of this 
region has been marked by a broken curve in these figures. The 
data beyond these boundaries is not believed to be accurate. 

The limitations noted above are not due to any failure of the 
Perturbation-Iteration Method but is due to the over simplified 
numerical procedure adopted to solve the integral equations. By 
using Galerkin ' s method or the method of least squares along with 
additional terms in the expansions fcr the charge and current 
distributions results that are valid over a wider range of 
parameter values could be obtained. 

In the region where the results become inaccurate the calculated 
current and charge distributions no longer have the expected 
behavior. This Is seen in Fig’s. 2.16 - 2.19 where graphs of the 
charge and current are displayed for sufficiently accurate and inac- 
curate expansions. The graphs only indicate one-half of the micro- 








strip as the entitles are symmetric around the origin. The Per- 
turbation- Iteration Method is generally valid for values of K e £f 
up to approximately .9*. The range of frequencies where K e f £ 

is valid also contains valid values of Z . 
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TABLE 



98.8 I 98.8 99.0 1 99.6 1101.9 1104.9 
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TABLE 11 


















































TABLE 13 



15.5 1 15.7 16.2 18.1 1 21.8 28.2 40.3 65.8 
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TABLE 17 



26.4 26.6 26.9 28.1 32.0 38.4 45.5 63.8 77.3 






































































TABLE 21 



113.0 















































TABLE 23 


























































TABLE 25 


















































TABLE 27 













































TABLE 29 





















































































TABLE 33 
































































TABLE 37 



60.2 60.3 60.6 61.6 64.5 67.3 69.9 72.5 75.6 
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Figure 2.4 Characteristic impedance Z c vs. frequency for tc = 2. — ratios 

indicated on graph. H = lmm. 




Frequency (GHz ) 











Figure 2.8 Characteristic impedance vs. frequency for k = 6. — ratios 

H 

indicated on graph, h = 1mm. 
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Figure 2.10 Characteristic impedance vs. frequency for ic - 8. — ratios 

H 

indicated on graph, h = 1mm. 
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Figure 2.14 Characteristic impedance Z vs. frequency for * = 12. jj ratios 
indicated on graph, h - Iran. 
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CHAPTER 4 


CONCLUSIONS 

This Chesls has developed a Perturbation-Iteration Method for 
analyzing microwave striplines. The theory began with the static 
analysis being developed followed by the higher order analysis. 

The static analysis was developed by using two static Green's 
functions. The higher order analysis was developed using three 
frequency dependent Green's functions and satisfying the frequency 
dependent wave equations. Extensive tables and graphs of computed 
values of and Z c over a range of frequencies was presented 

for different structures. Data was also presented which compared 
the method with previously developed methods and experimental data. 

The comparison showed that the results were accurate to at least 
1% over most of the parameter range of practical interest. 

The Perturbation-Iteration Method is a relatively easy method 
to use as it contains no transcentental dispersion equations that 
must be solved. Rather, it employs a set of equations which, when 
solved simultaneously, produce answers for 0 and Z c directly and 
which satisfy the field requirements of the structure. The ease of 
using this method is also enhanced by the fact that it calculates 
values based on ratios of the strlpllne dimensions. Therefore, 
once one set of dimensions are used for computations, other structures 
with the same dimension ratios can be directly analyzed by employing 
a frequency scaling. An example of this feature is a structure with 
mlcrostrlp width of 2 units, height above the ground plane of 1 unit 
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and frequency of 1 GHz has the same characteristics as a structure 
of mlcro8trlp width 1 unit, height above the ground plane of 0.5 
units, and frequency 2 GHz. This scaling feature enhances the 
value of the method greatly. 

Over the course of developing the method, another method was 

developed which proved inadequate. It used a perturbation series 

which expanded 0 in odd powers of frequency. It initially solved 

for a higher order solution. The static solution was identical as 

used in the Perturbation-Iteration Method. The higher order 

solution used the same static Green's functions as the static case 

along with frequency dependent sources. The method proved to be 

unsuitable as the computed values of 8 and Z did not converge 

c 

to any exact values because they were dependent on the dimension a. 
This dependence led to the abandonment of the metnod in favor of 
the Perturbation-Iteration Method. 

The validity of the numerical procedure used, holds over a 
broad range of frequencies but, becomes less accurate as the fre- 
quency is raised to high values. As the width of the microstrip 
is increased, the range of frequencies where the method is accurate 
decreases. However, the range of frequencies where the method 
hecomes less reliable is also the same range where higher order 
modes propagate down the strlpllne. These ranges of frequencies 
are not actually practical for use and thus are not a serious 
limitation to the effectiveness of the procedure for practical 
stripline structures. 
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The Perturbation-Iteration Method is a very practical method 

and can be applied in related areas. The theory can be modified 

to accomodate different structures such as slotlines and coplanar 

waveguides. The accuracy can be enhanced by adding extra cief- 

flcients to the charge and current density expansions which will 

make the theory valid over a broader range of frequencies. The 

Theory could be extended to separate and then solve for each mode 

of frequencies where multiple modes propagate. The method could 

also be reformulated to solve for 3 and Z in the Fourier 

c 

transform domain in a similar fashion as the fullwave analysis 
referred to in Chap. 1. It is believed that the loss in accuracy 
at very high frequencies is due to the use of a four term expansion 
of the charge nd current distribution on the strip along with 
point matching at four points along the strip. Other numerical 
procedures can be used to remove this limitation. Never the less, 
the simple procedure used is adequate for striplines of the type 
normally used in practice. 
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